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Sir: 

I, J. Nelson Wright, declare and state that: 

1 . I am a joint inventor, along with Lawrence J. Newell of the invention described 
and claimed in U.S. Patent Application No. 10/750,456 (the "Present 
Application"), filed on December 31, 2003. 

2. This Declaration establishes invention in this country of the subject matter in 
the currently pending claims in the Present Application before the earliest 
priority date of U.S. Patent No. 7,158,754 filed on July 1 , 2003 by Anderson 
("Anderson"). 

3. All of the work described within this Declaration was performed in the United 
States, by me or on my behalf and under my direction. 

4. I have reviewed my invention disclosure records, including the Exhibit 
submitted herewith, and readily conclude that the methods and systems as 
claimed in the Present Application were conceived and reduced to practice 
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prior to July 1, 2003, the filing date of the Anderson Patent. 

5. In support of this conclusion, I have attached Exhibit A, which represents a 
redacted invention disclosure document entitled "Array Calibration" 
("Invention Disclosure"). Although the dates have been removed from the 
enclosed Invention Disclosure, I represent that this document represents 
conception of the invention at least prior to July 1, 2003. In particular, a 
comparison between the figures of the Present Application and the figures of 
the Invention Disclosure illustrate that the figures in the Invention Disclosure 
correspond closely to Figures 4-9 of the Present Application. 

6. Further in support of this conclusion, I have attached Exhibit B, which 
represents a redacted schematic drawing of a working prototype of the array 
disclosed and claimed in the pending US Application 10/750,456 dated at 
least prior to July 1 , 2003. This schematic drawing is evidence of actual 
reduction to practice prior to the July 1 , 2003 critical date of Anderson. 

7. I further declare that all statements herein made of my own knowledge are 
true and that all statements made on information or belief are believed to be 
true; and further, that the statements are made with the knowledge that the 
making of willful or false statements or the like is punishable by fine or 
imprisonment, or both, under Section 1001 of Title 18 of the United States 
Code, and may jeopardize the validity of any patent issuing from the present 
patent application. 
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Array Calibration 

Objective: 

The objectives of this writeup are to: 

1) Reformulate Nelson's writeup of calibration from 1 1/19/02 in terms of 
differential and common-mode calibration operations. The primary 
motivation for doing this is to reduce the number of calibration operations 
which must be executed dynamically (up to a four-fold reduction). 

2) Extend Nelson's writeup to include treatment of additional pertinent parasitic 
elements or error sources to be incorporated for accurate calibration. 

3) Specify the calibration algorithm and EEPROM contents in sufficient detail to 
allow software developers to code the algorithm. 

Outline: 

The outline of this description is follows: 

1) Calibration Architecture 

2) Single-ended vs. Differential-Mode/Common-Mode Representation 

3) Model of Array Admittances 

4) Dynamic vs. Static Calibration 

5) Mitigation of channel-channel mismatch in voltage calibration sources 

6) Effects of Differential Imbalance in the voltage calibration sources 

7) Effect of feedthrough and crosstalk in the Calibration Drive Circuitry 

8) Parameters to be stored in EEPROM. 

9) Summary of the calibration algorithm. 

1) Calibration Architecture: 

Figure la is a schematic of a single differential preamplifier and associated calibration 
sources for the Pharos sense array. The preamplifier and array are discussed in detail 
elsewhere. The precision calibration sources are implemented by switching square-wave 
waveforms into any of four injection ports per differential preamplifier channel. The 
two injection ports through large-valued, precision resistors Rref make up the Ical, or 
current injection calibration ports. The two injection ports through precision resistor 
dividers (Rdivi and Rdivj) constitue the Vcal, or voltage calibration ports. The square- 
wave waveforms are generated with registers implemented in the AC family of CMOS 
logic. This family of logic switches the register outputs to either GND or V s through 
very low-valued switch resistances. A stable and accurate supply voltage Vs guarantees 
that square-waves of precise amplitude and phase are generated. The calibration sources 
can be activated in any desired pattern, either singly, differentially, or in common-mode. 
The receiver board contains filtering to eliminate all higher harmonics leaving only the 
fundamental of the calibration waveform. All signal-path circuitry is operated in a highly 
linear manner. 
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Figure 1a Figure 1b 



Figure lb is an equivalent circuit of the single differential preamplifier from figure la 
which can be derived using the hyprid-pi model assuming no impedance between the 
collector and base terminals of the transistor. The square-wave calibration sources have 
been replaced by idealized sinusoidal voltage and current sources with rms amplitudes 
related to the supply voltage V s as shown in the figure. 

Figure 2, depicts the entire M=32 coil array, N=64 current calibration sources, N voltage 
calibration sources and the network consisting of N/2=32 differential preamplifiers. As 
indicated by the dashed lines in figure lb, for analysis purposes, the Rref resistor, 
snubber and biasing network are incorporated in the Y L array matrix. The Vext vector 
represents the induced open-circuit Faraday voltage in the M sense coils. The NxM Y T 
'transadmittance' matrix relates the N short-circuit currents II measured at the calibration 
current injection points to the M elements of Vext- NxN Yl represents the 'self- 
admittance' of the array as measured from the calibration current injection points. Note 
that the N=2*M 2-ports on the left of the array can be considered as single-ended signals 
referred to a common reference, or as M differential ports and M common-mode ports 
referred to a common reference (with consistent representations of Yl and Y T ). 
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Preamplifiers, cables and 
Receiver board analog 



V =7 *l 

v p i- p i p 

V 0 =Z t *'p 



is* 



Coils, Pi 
circuitry, Parasitics 



" L =Y L *V L (Vext = 0 ) 
l L =Y T *V EXT ( v l= o, : 
OR 



Similarly, the preamps, cables and receiver board are characterized by an NxN (64x64) 
complex 'transimpedance' Z T matrix, an NxN complex 'self-impedance' Z P matrix, Nxl 
output voltage vector V 0 , Nxl input voltage vector V P and Nxl input current I P . 

The N calibration current sources Ical and N calibration voltage sources Vcal are related 
to the array and preamp voltages and currents by the nodal equations shown in the figure. 

In normal (localization) operation, V C al=0 and I C al=0 and the output voltage vector Vo 
and the 'hypothetical' short-circuit current II sc can be found using the equations in the 
figure from the Vext voltages: 

= • (I MxM + Y L -Z P ) • Y T • 

I =Y -V 

J L_SC 1 T EXT 

= -T-V Q 

Thus multiplying the measured output voltage vector by T gives the short-circuit current 
measurement of the array. T 1 is defined as follows and is measured during calibration: 



T- } -z T .(i MxM+ Y L .z P r 



During excitation of the calibration currents, with Vcal =0 and Vext =0: 
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*0 - %T ' (I*hM + Zp) ' " ^ CAL 

= T~'-I CAL 

During excitation of the calibration voltage, with Ical =0 and V EX t =0: 

Yo=-Z T -(I MM +Y L -Z p r.Y L -V CAL 
= -T- 1 -Y L -V CAL 

A calibration sequence consisting of measuring column vector V 0 while stepping through 
each of the Ical elements one at a time measures the columns of T 1 . Stepping through 
each of the Vcal elements one at a time measures the columns of T l *Y L . This 
architecture then allows us to measure Yl, and to measure the short-circuit current 
resulting from Vext: Il_sc= Y t *Vext- All we have to do now is relate the measured 
NxN Y L to the needed NxM Y T , and we can determine, with high accuracy, the original 
Vext- 

2) Single-ended representation vs. Differential-Mode/Common-Mode 
representation: 

As we have a fully differential circuit topology throughout the array and receive signal 
chain, there are advantages to representing the Z, Y and T matrices noted above with 
M=32 differential and M=32 common mode ports, rather than N=2M=64 single-ended 
ports. Each port always consists of two terminals. The single-ended and common-mode 
ports always use the common ground reference as one terminal, the differential ports do 
not use a ground reference. 

Denoting I+, L, V+ and V. as 32x1 column vectors corresponding to the positive and 
negative single-ended terminal currents and voltages for the 32 pairs of ports of our Y 
matrix, we can create the 64x1 column vectors and 64x64 Y matrix segmented as 
follows: 

Alternately, we may perform the coordinate transformation to differential and common- 
mode currents and voltages: 

' D =(W-)/2 
/ c =(/ + +/.)/2 

v D =v + -v_ 
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VI = \ y dd \ Yoc jv*] _ Lr ik±2kzjkzL± 1 jkzjk±jVzjk 1fc1 



Y DD represents the short-circuit differential current for a differential voltage drive. 

Ycc represets the short-circuit common-mode current for a common-mode voltage drive. 

Ydc and Ycd represent the mode coupling admittances. 

3) Model of Array Admittances: 

This section will discuss the relation of the measurable Y L (admittance measured from 
the array coil terminals) to the desired Y T (the transadmittance from the Faraday induced 
voltage to the array coil terminals). 

A number of different intentional circuit elements and unintentional (parasitic) circuit 
elements contribute to the total Y L . It can be shown that the admittances from these other 
elements just add to the desired transadmittance, Yj. Determining Yt from the measured 
Y L just involves determining and subtracting off the undesired components of Y L . 

Three types of array admittances can be identified with respect to calculating the desired 
Y T . 



3.1) The 'ideal' sense coils themselves: These components have no common mode 
admittance and no common to differential or differential to common mode 
coupling terms. The differential admittance is just the M.* 1 matrix described in 
Nelson's writeup. Note that the 64x32 Y T is just the first 32 columns of Y L i . 



3.2) Circuit elements which are connected directly between array coil terminals and 
the common reference ground potential: These components have a differential- 
mode admittance that is equal to the common-mode admittance. The differential 
mode components can be determined by measurement of the common-mode 
components of Y L . Examples in this category include: 

- Preamp bias resistors 

- Impedance of the Norton equivalent resistor for the current cal. source Rref 

- Stray capacitance to GND of differential pairs and coil 

All these components show up only on the diagonal of each block. In general, 
these components are reasonably well matched (-1%) and fairly small to start 



Y L =Y U +Y U+ Y L 
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with. If the shunt impedances to GND are denoted z lp , z, n , z 2p , z 2 „, ... z Np , z N „, 
then the admittance matrix looks like: 



0 



0 



0 



0 



0 



0 



3.3) Those elements whose common-mode and differential-mode admittances are not 
the same, and hence must be characterized in engineering or in the factory and 
included in the EEPROM. There are two types here, which will be modeled 
separately. 



The first type, Y UA covers the differential shunt impedance across each coil. This is 
from the intentional snubber network, and the unintentional parasitic capacitance across 
the coil and differential pairs leading from the coil. If these individual impedances are 
denoted z 3J , z 32 , ... z 3 n, the admittance matrix looks like: 



0 0 
0 0 



0 0 



0 0 
0 0 



0 0 



0 0 
0 0 



0 0 



The second type, Yl3b , covers the small capacitive coupling between the differential 
pairs and coils in the array. Unfortunately, this is not expected to appear mainly 
common-mode or differential mode. However it will have O's on the diagonals of each 
block. This type will need to be characterized in differential mode and stored in 
EEPROM. 
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Thus to determine the desired Ytdd^Mc" 1 we take the measured differential Y LD d, 
subtract the Yuadd and Yubdd derived from the EEPROM contents, and subtract the 
diagonal elements of the measured common-mode YLccaiagonai, 

Note that only the diagonal elements of the measured Ylcc should be subtracted off. The 
off-diagonal elements include common-mode capacitance (Yubcc) between the coils that 
do not affect Yt- 

4) Dynamic vs. Static Calibration: 

Some of the array and receiver parameters discussed above are subject to variation with 
temperature (such as the preamp gain and input impedance represented by Z T , Zp and T). 
Other parameters such as Yt, Y l or Mc are subject to variation on much shorter time 
scales as metal objects may be introduced to or removed from the vicinity of the array. 
Still others are expected to remain essentially static over time, such as Ylcc- It is 
therefore important to execute measurements of the time-varying parameters interleaved 
with localization measurements to assure accurate calibration of the Pharos system. 

Executing calibration measurements frequently, however, reduces the percentage of time 
during which localization measurements can be made, which in turn will impact 
adversely either the update rate or variance of the localization measurements. 

As only one block of the Yl matrix, Yld D , is prone to variation, we can save considerable 
time by exciting the VCAL sources in differential mode only, and making only the 
differential receive measurement. This reduces by 75% the time required to make the 
voltage calibration measurements. 

A similar 75% savings can be achieved in the current calibration measurements if the T 
matrix is block-diagonal. In the calibration process outlined in section 1, the columns of 
T 1 are measured sequentially, and the desired T matrix is calculated with a matrix 
inversion. In general, the full T l matrix must be measured to determine T, but if T 1 is 
block-diagonal, the diagonal blocks can be measured and inverted independently. 

Tinv DD 0 T'Trmv-i, 0 1 
0 Tinv cc \ [ 0 71mg.\ 

In the Pharos system, measurements of the T 1 matrix indicate that the off-diagonal 
blocks are 40 to 50 dB smaller than the diagonal blocks. The error introduced in making 
the block-diagonal assumption in the matrix inversion goes as the square of the relative 
magnitudes of the off-diagonal blocks. We can then expect errors of less than 100 ppm 
in making this assumption. 
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Thus only differential -mode excitation and differential-mode receive measurements need 
to be made dynamically. The full calibration need only be executed on a much less 
frequent basis, such as once per patient, once per day, or even just once in the factory. 

5) Mitigation of channel-to-channel mismatch in voltage calibration sources: 

In characterizing Y L one column at a time, the 'short-circuit' current is measured while 
driving one channel (differentially, or common-mode) at a time. Mismatches in the 
individual (single-ended) calibration voltage sources result in the following effects: 

• Channel-channel amplitude mismatch in the differential drive voltages 

• Channel-channel amplitude mismatch in the common-mode drive voltages. 

• Crosstalk into the common mode when driving differentially 

• Crosstalk into the differential mode when driving common-mode. 

Crosstalk between modes is treated in the next section. 

Mismatches in the current calibration sources are expected to be limited by the accuracy 
of the Rref resistor, or about 0.02%. Achieving this type of accuracy in the VCAL 
voltage dividers, though, is not feasible. 

Mismatches in the voltage divider unique to each channel will result in gain errors 
common to each column of the measured admittance. Let Y M be the measured 
(corrupted) admittance, and Gv* 1 be the (complex) diagonal matrix representing the 
(unknown) gain errors that caused the corruption: 



We want to find that Gv = diag{[gi g2 g3 • ■ ■ g32]} to correct the measured Ym and achieve 
the true Yl. Reciprocity will guarantee that Yl is a symmetric matrix. A least squares fit 
can be applied to 'back out' the differential voltage drive mismatches that were present. 
The equations describing reciprocity of Yl are: 



This is a system of V£*(32 A 2-32)=496 equations with only 32 (complex) unknowns (g] 
... g32). Of course the least squares fit to this system of equations will be gi=g2= ... g32 = 
0. To assure the solution we want, we need to add one more equation that assures the 
mean of the gains is 1 : gi+g2+ . . . grc = 32. 




Z.31 



L2I 



or 



Y MnSl ~ Y M2)&) 
Y MnSi ~ Y Ml\%\ 



This system of equations can be captured in matrix form as: 
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with: 



Ag = b 

g = kx g 2 & Ai gitl 



A = [0 0 0 0 0 0 0|32j 
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The least squares solution to this overdetermined system of equations is: 

Any differential-mode channel-channel mismatch should also be present in the common- 
mode measurement. The same voltage source error fit obtained to achieve a symmetrical 
Y L in the differential block should be applied to the common-mode block as well. 



6) Effects of Differential Imbalance in the calibration drive voltage sources Vcal: 

Crosstalk from the differential mode to the common-mode can be analyzed as follows. 
When driving channel 1 differentially: 

i -_r v -v v - Xm I jog 1 Vcal_d [ 

h_sc- 1 Vo-'l v CA L -\y r~p \~y 

L CD I 'CCj \_ r CAL_CJ 

VcAi^fy 0 - 0 1 0 ■•• Oj 



instead of measuring just the first column of Ydd, we will also get e times the first 
column of Ydc added to it. 

Y D c contains only the mismatch in Yu shunt components, and Yub stray capacitances. 
For a rough order of magnitude estimate ... the shunt elements in Y L 2 have 1/10 the 
admittance of the diagonal elements of Ydd- The mismatch in these elements (which is 
what Ydc contains) is probably < 1/30 times that. Finally, e is on the order of 1/1000. 
This makes this error term on the order of 3 ppm of the diagonal elements of Y D d, and 
should be negligible. 
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Crosstalk from common-mode to differential mode will be the same order of magnitude. 

7) Effect of Feedthrough and Crosstalk in the Calibration Drive Circuitry: 

There are 4 types of feedthrough from calibration drive circuitry to the Vo output to 
consider: 

• Parasitic paths coupling the current-drive waveform into the V 0 output - diagonal 
terms. 

The dominant path expected here is the small amount of capacitance across the 
Rref resistor. Measurements indicate this is below 0. 1 pF. At 500 kHz, this 
results in a gain error below 1 part per 10,000. Also, the layout of each channel is 
identical, so this gain error should be the same across all channels. This effect 
will be ignored (but verified in the hardware). 

• Parasitic paths coupling the current-drive waveform into the Vo output - off- 
diagonal terms. 

Every attempt was made to minimize the possibility of this. For now, we will 
assume we don't have to compensate for this effect. This will be verified on the 
hardware by just removing the Rref resistor and measuring the S-vectors received 
in current calibration mode. 

• Parasitic paths coupling the voltage-drive waveform into the Vo output - diagonal 
terms. 

A significant effect is expected here from the collector-base capacitance in the 
preamp PNP transistor. We have not yet attempted to incorporate this in the 
circuit model, but errors from this source should be effectively mitigated by 
measuring and subtracting off the common-mode array admittance. In measuring 
the common-mode array admittance, the collector-base capacitance will appear as 
a 'negative' capacitance to GND in the array. This effect is present to the same 
degree in both differential and common-mode drive, so subtracting YuxDiagonai 
from the measured Yldd should be effective in compensating for this. 

• Parasitic paths coupling the voltage-drive waveform into the V 0 output - off- 
diagonal terms. 

As with the off-diagonal current-drive terms, we're hoping for the best on this one 
for now. 



There is also the potential for crosstalk between current-drive and voltage drive sources 
or between channels. Again, every attempt has been made to minimize this in the layout, 
and we will not attempt to treat this at this time. 
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8) Summary of EEPROM stored data: 

EEPROM Stored Data: 

Rref: Current calibration reference impedance 

Single real 16 bit value (in ohms) 

Nominal Value = 49900 (49.9 Kohms) 
Kdiv: Differential Voltage calibration divisor ratio. 

Single real 16 bit value (xle-6) 

Nominal Value = 6350 (0.006350) 
R1JI3: Sense Coil snubber/bias resistance values. 

Two real 16 bit values (in ohms) 

Nominal Values = 1124,5340 
Csnub: Sense Coil snubber resistance value 

Single real 16 bit value (in 0.01 pF) 

Nominal Value = 3900 (39 pF) 
C3A: Parasitic differential shunt capacitance across sense coils 

32 real 16 bit values (in 0.01 pF) 

Nominal Value = 500 (5 pF) 
C3B: Differential mode parasitic capacitance between sense coils 

-100 non-zero values of upper-triangular 32x32 real matrix of 16 bit 
values (in 0.01 pF) (nearest neighbors only) 

Nominal Value = 500 (5 pF) for nearest neighbors 
YLDDdiagronai: Results of dynamic cal in the factory at 300kHz, 400 kHz, 500 kHz 

32x3 complex 16 bit values (in units of 100 nS) 

Nominal value at 300 kHz = 6835 - j*47240 (6.835e-4 - j*4.724e-3 S) 

Nominal value at 400 kHz = 3879 - j*35751 (3.879e-4 - j*3.575e-3 S) 

Nominal value at 500 kHz = 2493 - j*28721 (2.493e-4 - j*2.872e-3 S) 
YLccdiagronai: Results of static cal in the factory at 300kHz, 400 kHz, 500 kHz 

32x3 complex 16 bit values (in units of 100 nS) 

Nominal value at 300 kHz = 874 - j*377 (8.74e-5 - j*3.77e-5 S) 

Nominal value at 400 kHz = 874 - j*503 (8.74e-5 - j*5.03e-5 S) 

Nominal value at 500 kHz = 874 - j*628 (8.74e-5 - j*6.28e-5 S) 

Calculation of Yua and Yun from EEPROM Stored Data: 

Let s=j*2*pi*Fc where Fc is the operating center frequency (unique for each beacon). 

For Yua, first we calculate the differential admittance of the snubber/bias network, and 
subtract off the common mode admittance of the snubber/bias network (which is 
accounted for in the common-mode measurement). This is a single parameter for all 
channels. 
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To get Yl3a, we then add the per channel capacitances, C 3A . 

Y L3A = Y SNUB ^32*32 + ^g{sC 3A ) 

Yub is just s times a capacitance matrix. This matrix is calculated by populating the non- 
zero elements and adding its transpose. (The pattern of non-zero elements to be updated 
later). 

9) Summary of the Calibration Algorithm: 

All calibration operations discussed are divided into two software processes: Static 
Calibration and Dynamic Calibration. Static calibration is currently planned to be 
executed at 'Session test' at the beginning of treatment for each patient. Dynamic 
Calibration is interleaved with localization measurements. The figures below indicate the 
processes involved with each. The individual blocks are described in more detail 
subsequently. 
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The next block describes how dynamic calibration measurements are interleaved with 
localization measurements. Note that N will be typically around 4. Each localization 
measurement will require roughly 100 msecs. Each calibration block will require 
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roughly 40 msec. Executing the full dynamic calibration cycle will then take roughly 4 
seconds. 



Interleaving of Dynamic Calibration Measurements and 
Loclalization Measurements 
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Receiver Gain (Gn) Calibration 

The objective of Receiver Gain Cal is to back out gain errors associated with the higher 
gain typically used in localization (Gr' 1 ). Calibration is performed at the low-gain 
setting. If localization is performed with the low-gain setting, no correction is needed. If 
localization is performed with the high-gain setting, the measured S-vectors should be 
left-multiplied by Gr prior to any other processing. 

The following steps are performed in receive gain cal: 

1) Configure receiver front-end for high-gain with inputs selected to V CM 

2) Measure S_HIGH complex values with cal-kernel, NMSI=1 
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3) Configure receiver front-end for low-gain with inputs selected to V_CM 

4) Measure S_LOW complex values with cal-kernel, NMSI=1 

5) Form grjnv, a 32x1 complex vector by taking the element-by element ratio of 
S_HIGH over S_LOW. 

6) Calculate Gr = (diag{gr_inv}) _1 

7) Output Gr 

Calibration Voltage Mismatch (Gy) Calibration: 

The objective of this process is to determine the channel-channel gain mismatches 
present in the calibration voltage sources. 

The following steps are performed to determine Gv: 

1 . Configure receiver in low gain and differential reception 

2. Excite Cal Current sources one at a time differentially. The 32 sets of 
measured S vectors form the matrix S_ICALdd. 

3. Excite Cal Voltage sources one at a time differentially. The 32 sets of 
measured S vectors form the matrix S_VCAL D d. 

4. Calculate T D d> YLDDMeaswed, as follows: 



DD ~ VlCAL DD J 

REF 

^LDDMeasured ~ ~~Z ^DD ' ^VCAL DP 
DIV 

5. Form the matrix A from the elements of Y L DDMeasured and construct the vector 
b as described in section 4. 

6. Obtain the correction gain vector g and matrix Gv as follows: 

g = ^A) l A"b 
G y =diagfe} 

7. Output Gv. 
Common-Mode Calibration: 

The objective of this process is to measure the diagonal terms of the common-mode 
block of Y L , designated as Y L ccdiagonai. 

The following steps are performed to determine Y L ccDiagonai : 

1) Configure receiver in low gain and single-ended receive (+ terminal) 
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2) Excite Cal Current sources one at a time in common-mode. The 32 sets of 
measured S vectors form the matrix S_ICAL C+ . 

3) Repeat steps 1 and 2 for - terminal to collect S_ICAL C - 

4) Excite Cal Voltage sources one at a time in common mode receiving single- 
ended both + and - terminals. The two 32x32 sets of measured S vectors 
form the matrices S_VCAL C+ and S_VCAL C .. 

5) Calculate T C c and use Gv (see below) to calculate Y L ccDiagonai as follows: 

D1V 

Zero out all non-diagonal elements of Ylcc to form YLccDiagonai 

6) Output YLCCDiagonal- 

Differential-Mode Calibration 

The two objectives of differential-mode calibration are to: 

• Frequently update the differential block of the T matrix relating measured 
S values to short circuit current at the preamp inputs 

• Frequently update the transfer function from Vext to differential short 
circuit current at the preamp inputs (ie Mc" 1 )- 

The following steps are performed for Differential Cal: 

1 . Configure receiver in low gain and differential reception 

2. Excite Cal Current sources one at a time differentially. The 32 sets of measured S 
vectors form the matrix S ICALdd. 

3. Excite Cal Voltage sources one at a time differentially. The 32 sets of measured S 
vectors form the matrix S_VCALdd- 

4. Calculate T D d, YLD D Measured, and M c as follows: 

REF 

^LDDMeasured = ~jl ^DD ' ^VCAL DD 
K Diy 



^LDDCometed ~ ^ LDDMeasured ' ^LCCDiagmal ^ Li. 

M c ~ YwDCormcted "2^ LDDConected ) 



5. Output Tdd and M c . 



Application of Calibration Results to raw S vectors for Localization: 

^CORRECTED = ' ^DD ' Gr ' ^RAW 
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